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S^ , Abstract. A model- independent experimental signature for flavor oscillations in the 

f^ ' neutrino signal from the next Galactic supernova (SN) would be the observation 

^D ■ of Earth matter effects. We calculate the probability for observing a Galactic SN 

^_J I shadowed by the Earth as a function of the detector's geographic latitude. This 

f"^ ■ probability depends only mildly on details of the Galactic SN distribution. A location 

^O . at the North Pole would be optimal with a shadowing probability of about 60%, but 

a far-northern location such as Pyhasalmi in Finland, the proposed site for a large- 
ly ■ volume scintillator detector, is almost equivalent (58%). We also consider several 
SpH| pairs of detector locations and calculate the probability that only one of them is 
O ' shadowed, allowing a comparison between a shadowed and a direct signal. For the 

South Pole combined with Kamioka this probability is almost 75%, for the South Pole 
^ I combined with Pyhasalmi it is almost 90%. One particular scenario consists of a large- 

volume scintillator detector located in Pyhasalmi to measure the geo-neutrino flux in 
a continental location and another such detector in Hawaii to measure it in an oceanic 
^% ' location. The probability that only one of them is shadowed exceeds 50% whereas the 

j^ ■ probability that at least one is shadowed is about 80%. We provide an online tool to 

calculate different shadowing probabilities for the one- and two-detector cases. 
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1. Introduction 

Galactic core-collapse supernovae are rare, perhaps a few per century tl,,2j. However, 
the large number of existing or future neutrino detectors with a broad range of science 
goals almost guarantees continuous exposure for several decades, so that a high-statistics 
supernova (SN) neutrino signal may eventually be observed. Such a measurement 
would provide a plethora of new insights that are crucial both for our astrophysical 
understanding of the core-collapse phenomenon [Hj and for using SNe as particle-physics 
laboratories j3]. Of particular interest would be the observation of signatures for flavor 
oscillations because this could address a key question of neutrino physics, i.e. whether 
the neutrino masses are ordered in a normal or inverted hierarchy IHl 13 IHl El QH] • 

The collapsed core of a SN emits neutrinos and anti-neutrinos of all flavors with 
comparable fluxes and spectra [HI [T21 HBl IHj- One expects differences between the 
spectra and fluxes of Ue and z/^^^ or between i/g and u^^r that are large enough to observe 
flavor oscillations, but it will be difficult to establish such effects solely on the basis of 
a z/e or z/e "spectral hardening" relative to theoretical expectations. Therefore, in the 
recent literature the importance of model-independent signatures has been emphasized, 
e.g. in association with the prompt z/g neutronization burst ^H] or with shock-wave 
propagation [121 HZl dl CH]- One unequivocal signature would be the observation of 
Earth matter effects, because they induce a characteristic energy-dependent modulation 
on the measured flux [HI El El 1201 1211 • ^^ ^ large- volume liquid scintillator detector 
such as the proposed 50 kt Low Energy Neutrino Astronomy (LENA) project |22] this 
modulation could be detected with high statistical signiflcance ^Hl- Moreover, even 
if single detectors can not resolve these modulations, comparing the signals from a 
shadowed with an unshadowed detector may allow one to diagnose Earth effects |Sj. 

These intriguing possibilities have motivated us to investigate the probability for 
given detector locations to observe the next Galactic SN in an Earth-shadowed position. 
Most of the Milky Way is in the southern sky so that a northern location is obviously 
preferred, but a quantitative determination of the probability distribution of Earth- 
shadowing as a function of geographic latitude is missing, except for a brief discussion in 
Ref . in] . An additional motivation for our work is that large- volume scintillator detectors 
are being discussed for the purpose of geo-neutrino observations. After KamLAND's 
pioneering measurement of the z/g flux from uranium and thorium /5-decays in the 
Earth |2S1 , performed in a complex geophysical environment with a large Ue background 
from power reactors, the exploration of different detector sites has become crucial. 
The choice of location could be influenced by the role of such detectors as excellent 
SN neutrino observatories, in particular because they need to operate for a long time 
to accumulate meaningful statistics for geo-neutrino studies. For example, if a large- 
volume scintillator detector were built in the Pyhasalmi mine (Finland) J21I to measure 
the geo-neutrino flux mainly from the continental crust and another one in Hawaii |2H] 
to measure neutrinos from the oceanic crust, the "geographic complementarity" of these 
locations with regard to SN shadowing is important. 
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We begin in Sec. |21 with a discussion of the distribution of core-collapse SNe in 
the Galaxy and then determine the SN probability distribution in the sky. In Sec. El 
we determine the probability for Earth shadowing as a function of geographic location 
of one or two detectors. We also provide an online tool for the calculation of these 
probabilities J2E]- We conclude our work in Sec. IH 

2. Supernova distribution in the Milky Way 

In this section we characterize the SN probability distribution in our Galaxy. In Sec. 12.11 
we present the models adopted for the SN volume distribution. In Sec. 12. 21 we discuss the 
probability of a core-collapse SN in terms of Galactic and Equatorial coordinates. Finally 
in Sec. l2.3l we estimate the distance distribution of SN events. The SN distribution in the 
Galaxy is expressed in cylindrical galactocentric coordinates {r,z,6), where the origin 
corresponds to the Galactic center, r indicates the radial coordinate, 6 the azimuthal 
angle and z the height with respect to the Galactic plane. 

2.1. Models for the volume distribution 

The probability distribution of core-collapse SNe in the Milky Way is not well known. 
These SNe mark the final evolution of massive stars and thus must occur in regions of 
active star formation, i.e. in the Galactic spiral arms. As proxies for the core-collapse 
SN distribution one can use either observations of other Galaxies, or in our Galaxy 
the distribution of pulsars and SN remnants (SNRs), the distribution of molecular 
hydrogen (H2) and ionized hydrogen (HII) and the distribution of OB-star formation 
regions. (For a review see Ref. [2Zj and references there.) These observables are either 
directly connected with core-collapse events (SNRs, pulsars) or with young, massive 
star formation activity and the related emission of ultraviolet light (H2, HII, OB stars). 
All of these observables are consistent with a deficit of SNe in the inner Galaxy and a 
maximum of the probability at 3.0-5.5 kpc galactocentric distance. 

The star formation activity is not smoothly distributed in the Galaxy. Besides 
generally following the spiral arms, it can be concentrated in small regions or spike-like 
complexes like the Cygnus OB association. Small regions of high star-forming activity 
have also been found within 50 pc from the Galactic center |2HI that may contribute up 
to 1% of the total Galactic star formation rate, although this finding does not seem to 
contradict the overall picture of a reduced SN rate in the inner Galaxy. 

However, we are only interested in the SN distribution projected on the sky and the 
Earth's rotation introduces an additional averaging effect. Therefore, it will be enough 
to consider a smooth distribution with azimuthal symmetry. In particular, we shall use 
the following common parametrization for the Galactic surface density of core-collapse 
(cc) events, 

<7cc{r) ocr^exp(-r/u) , (1) 
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Figure 1. Galactic surface density of core-collapse SNe as a function of galactocentric 
radius according to Eq. (^. Solid curve: Our benchmark case with the parameters of 
Eq. ^. Dashed curve: Alternative parameters of Eq. Q. Dotted curve: Type la SNe 
according to Eq. ©. The surface densities are normalized as j cr{r) 2tt rdr = 1. 



(2) 



where r is the galactocentric radius. For the birth location of neutron stars, a fiducial 
distribution of this form was suggested with the parameters 

{ f = A 

Neutron stars: < ' 

[ n = 1.25 kpc . 

These parameters are consistent with several SN-related observables, even though large 
uncertainties remain jSHj- On the other hand, the pulsar distribution indicates jHOj 

Pulsars: U = 2-35 , 

\ u = 1.528 kpc. ^ ' 

We will use the parameters of Eq. Q as our benchmark values and those of Eq. (jH)) 
as an alternative model to illustrate the dependence of our results on different input 
choices. 

In Fig. C] we show the normalized Galactic surface density of core-collapse SNe 
according to Eq. (P) as a function of the galactocentric distance for the two choices of 
parameters given in Eqs. and (jH))- Our benchmark distribution shows a peak at 
r = 5 kpc, while for the parameters of Eq. Q the surface density peaks at a lower 
distance r ~ 3 kpc. The distributions are normalized as /drcr(r)27rr = 1, i.e. the 
surface densities a are given in SNe per kpc^. Of course, to obtain the number of 
SN events per Galactic unit surface and unit time, the quantities a would have to be 
multiplied with the integrated Galactic SN rate. However, since we are only interested 
in the relative probability of SNe in different regions of the sky, the overall rate is not 
important for our study. 
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The vertical distribution of neutron stars at birth with respect to the Galactic plane 
can be approximated by the superposition of a thin Gaussian disk with a scale height 
of 212 pc and a thick disk with three times this scale height, containing 55% and 45% 
of the pulsars, respectively p7| . 



7^cc(-2) oc 0.79 exp 



+ 0.21 exp 



636 pc 



(4) 



V212 pc^ 

where z is the height above the Galactic plane. We assume this vertical distribution to 
be independent of galactocentric distance j3Q| so that 

Uccir, z) oc acc(r) TZcc{,z) (5) 

is the volume distribution. 

We stress that the distribution of core-collapse SNe differs significantly from the 
overall matter distribution, particularly in the inner part of the Galaxy. The distribution 
of Type la SNe — that are believed to originate from old stars in binary systems — more 
closely follows the matter distribution. It can be parameterized as |2Z| 

nia(r, z) = a,.{r)nUz) oc exp (-^;^) exp (-^) • (6) 

In Fig. C] we show the monotonically falling SNe la surface density aia.{r) for comparison 
with the core-collapse case. 

2.2. Projection on the sky 

The probability of a core-collapse SN as a function of the Galactic longitude / and 
latitude b is given by an integration along the line of sight, 

P{l,b) (X dsncc[r{s,l,b),z{s,b)] , (7) 

Jo 



/OO o \ / 

r = is cos b + drr, — 2 s drr, cos / cos b 



where 

z = s sinb . (8) 

Here, — vr < I < it, — 7r/2 < b < tt/2, and rf© ~ 8.5 kpc is the solar distance from the 
Galactic center. The function P{1, b) can be recast in terms of the Equatorial coordinates 
(a, (5), in which the reference plane is the Earth equator instead of the Galactic plane. 
The transformation relating the two systems of coordinates is 

sin 5 = sin 6 sin 5ngp + cos b cos 5ngp sin(/ — /q) , (9) 

cos(a — ao) = cos(/ — /o) cos6/cos5 , (10) 

sin(a — ao) = [— sinfocos^NCP + cos6sin(5NGpsiii(^ — /o)]/cos5 , (11) 

where, for the Julian epoch J2000, the coordinates of the north Galactic pole (NGP) are 
ttNGP = 12 h 51.42 m and (5ngp = 27°07.8'. Hence, the ascending node of the Galactic 
equator is at Oq = 282.86° and /q = 32.93°. Further details on astronomical coordinate 
systems can be found, for example, in Ref. p?T| . 
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Figure 2. Probability distribution P{a, S) of core-collapse SNe in the sky. Top: Earth 
Equatorial coordinates. Bottom: Average over right-ascension, see Eq. H12|l . 



The sky map in Equatorial coordinates is shown in the upper panel of Fig. |21for our 
benchmark distribution. However, since the neutrino detectors are fixed to the Earth, 
we only need this distribution averaged over time, i.e. over right ascension a. Therefore, 
all relevant information is contained in the "exposure probability function" 



/+'ir 
da P{a, 6) . 
-IT 



(12) 



It provides the probability distribution of the arrival direction of a SN signal in terms 
of declination. We show the time-averaged sky map in the bottom panel of Fig. |21 
normalized as J d6uj{6) cos6 = 1. This normalization also fixes the overall constant 
for the function P{a, 6) plotted in the top panel of Fig. |21 Fig. |21 clearly shows the 
preference of southern locations in the sky and also shows the polar regions that are 
nearly void of SNe because the Milky Way extends approximately between declinations 
of -60° to +60°. 
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Figure 3. The "exposure function" llj{5) for different Galactic SNe distributions. Solid 
curve: Bencfiniark model based on Eq. (Q with the parameters of Eq. (j^J. Dashed 
curve: Same with the alternative parameters of Eq. |(2Jl. Dotted curve: Distribution of 
SNe la for comparison. 



In Fig. iniwe show the exposure probabihty function 00(6) for different assumptions 
about the Galactic SN distribution, i.e. our benchmark distribution Eq. (Q) with the 
parameters of Eq. ((2)) , the alternative parameters of Eq. (jHl) , and as an extreme case the 
SN la distribution of Eq. © which, of course, is not realistic for core-collapse SNe. The 
normalization / dSuj{S) cos 6 = 1 is adopted throughout. In this normalization oiuj{S) an 
isotropic SN distribution would correspond to a horizontal line in Fig. El The exposure 
function depends only mildly on details of the assumed Galactic SN distribution, because 
the geometric effect dominates that the SN probability is negligible outside the Galactic 
disk. The largest model variation occurs around a declination of —30° near the Galactic 
center region. The edges at about ±60° correspond to the region where the Milky Way 
would end in the sky if it were infinitely thin. The tails beyond these declinations come 
from the vertical extension of the Galactic disk around us where we have assumed the 
solar system to be located exactly in the Galactic plane. 



2.3. Distance distribution 

As an aside, we use our assumed core-collapse SN distribution in the Galaxy to evaluate 
the SN distance distribution relative to the solar system. The average distance is 

/ ncc{r, z)d{r, z, 6) r drdzdO 



where 



(4c) 



d{r,z,9) 



J riccir, z)2'n'rdrdz 

{X - Xq)^ + (y - 2/0)^ + (Z- Zq)'^ 



(13) 



1/2 
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Figure 4. SN probability vs. distance from the Sun. Solid curve: Benchmark model 
based on Eqs. (Q) and |21l. Dotted curve: SNe la assuming Eq. ©. 



+ z + cIq — 2rdQ cos 6 



1/2 



(14) 

For the distribution rice of Eq. (0) with our benchmark parameters of Eq. (J2|) we find 
(dec) = 10.7 kpc with a rms dispersion of 4.9 kpc. While the average distance agrees 
with the fiducial distance of 10 kpc that is frequently assumed in the literature, we note 
that the dispersion of distances is very large. This is clearly visible in Fig. |3] where 
we plot the probability distribution Il{d) of SN events as a function of the distance 
d from the Sun, so that /^ Il{x)dx gives the probability that a SN happens between 
distance a and b from the Sun. The dip in the middle corresponds to the deficit of core- 
collapse SNe in the Galactic center region. For comparison we also show the probability 
distribution for SNe la according to Eq. (jHl). Their average distance is 11.9 kpc with a 
rms dispersion of 6.0 kpc. On the other hand, the median distances are quite similar for 
the two distributions, and equal to about 10.9 kpc. The large-distance behavior for both 
cases differs significantly due to the different scales in the exponential tails of Eqs. (Hj) 
and (jHI)- Of course, the scarcity of data at large galactocentric distances implies that 
the extrapolation may be unphysical. 



3. Optimal detector location 

3.1. One detector 

Armed with the exposure function uj{6) shown in Fig.lHlwe now determine the probability 
that a detector located at a geographic latitude A will observe the next Galactic SN 
below the Earth horizon. Of course, observable matter effects would require a minimal 
path-length of a few thousand kilometers ^^ so that our Earth-shadowing criterion is 
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somewhat schematic. Moreover, the phenomenological signature of Earth matter effects 
also depends on whether or not the neutrinos cross the core J2I1 E21 (HH] ■ We use a core 
radius -Re = 3486 km and an Earth radius of R^ = 6371 km j31]. The Earth or core 
shadowing condition for a source with altitude a with respect to the horizon is 

{0 , Earth shadowing , 

(15) 
— sin Oc , Core shadowing , 

where 



sin ac = ^1 - {RclR®Y = 0.837 . (16) 

In general, for a neutrino path-length L in the Earth, one has 

K = —smaL = • (17) 

The altitude of an object of equatorial coordinates (a, 5) is 

sin a = sin A sin 6 + cos A cos 6 cos H , (18) 

where H = t—a is the hour angle and t the local sidereal time. Since no time information 
on the next SN is available in advance, H has to be taken as a random variable. 

For polar locations (A = ± 7r/2) the shadowing condition simplifies, because 
cos A = and the time dependence disappears. The geometrical probability ^^(A,^) 
for a trajectory to satisfy one of the shadowing conditions is 

P.{±l,S^ =e(«:Tsin5), (19) 

where G(x) is the usual step function. A similar simplification holds for objects located 
at the celestial poles [6 = ± 7r/2) for which sin a = it sin A and 

p«(^A,±0 =0(ft:TsinA) . (20) 

Equation (j^ is related to Eq. ^T^ by the symmetry p^{X,6) = Pk{S,X). This is a 
general property following directly from Eq. (fTHj) . 

Apart from these special cases, we have both cos A > and cos 5 > so that the 
general shadowing condition Eq. ()15p can be recast as 

cos i7 < — tan A tan (5 H -, (21) 

cos A cos 

or equivalently 

cos if > tan A tan 5 r - = hJX,5). (22) 

cos A cos 

Note that we have replaced H + 7t with H in the argument of the cosine, since both are 
random variables. Unless Eq. ()22j) is always or never satisfied, one has 

— arccos/iK(A, (5) < H{m.od 27r) < aiccos hf^{X, 6) . (23) 

The general solution is then 

0[1 - K{\, 5)] e[-K{X, 5) - 1] for \K{\, 5)\>1, 

P'^(>^^^) = { 1 , ,. ,, , . (24) 

— arccos /Ik(A, o) otherwise. 
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Figure 5. Probability for Earth sliadowing (upper curves) and core shadowing (lower 
curves). Solid lines: Galactic benchmark distribution with parameters given by Eq. |2Il. 
Dashed lines: Alternative parameters of Eq. Q ■ The latitude of the Pyhasalmi site is 
indicated by the vertical dotted line. 



This follows because all values of the random variable i7(mod27r) in the interval (— tt, tt] 
are equally likely. As a check of Eq. (^^ we note that for a detector at the equator 
(A = 0) it simplifies to 



Pk(0, 6) = — arccos 

n \ cos 



K 



(25) 



reproducing the trivial result po{0,6) = 1/2. 

In order to obtain the Earth shadowing probability for the neutrinos from the next 
Galactic SN, given the detector latitude A, we must convolve with the distribution of 
source declination angles, i.e. with the exposure function, so that 

P.(A) = Jd6 cos 6 p.(A, 6) uj{6) , (26) 

where we have used the normalization J d6 cos6uj{6) = 1. 

We show the Earth and core shadowing probability as a function of detector latitude 
in Fig. 13 The solid lines refer to our Galactic benchmark distribution Eq. (Q) with 
parameters Eq. (j21) whereas the dashed lines refer to the alternative parameters Eq. (jS)) . 
Once more we find that the dependence on details of the Galactic distribution is mild. 
Note that for the Earth-shadowing case -Po(A) + Po{—X) = 1 because all information on 
the longitude is lost and two locations at A and —A are complementary: if one detector 
is shadowed, the antipodal one is certain not to be shadowed. 

The largest probability for Earth shadowing is a the North pole, but the Pyhasalmi 
site in Finland is almost equivalent. The behavior for the core shadowing condition is 
similar, although the advantage of a northern site is somewhat boosted: for k = — sin Oc, 
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Table 1. Representative locations of proposed or existing SN neutrino detectors and 
neutrino shadowing probabilities, assuming our benchmark Galactic distribution. 



Location 


Latitude 


Longitude 


Shadowing 


probability 












Earth 


Core 


Pyhasalmi, Finland 


63.66° 


N 


26.04° 


E 


0.581 


0.116 


Soudan, USA 


47.82° 


N 


92.23° 


W 


0.572 


0.112 


Frejus, France- Italy 


43.43° 


N 


6.73° 


E 


0.568 


0.110 


Kamioka, Japan 


36.27° 


N 


137.3° 


E 


0.560 


0.104 


Hawaii, USA 


19.70° 


N 


156.30° 


W 


0.528 


0.082 


Sydney, Australia 


33.87° 


S 


151.22° 


E 


0.445 


0.069 


South Pole 


90° 


S 


— 




0.414 


0.064 



the probability P^i^) varies by almost a factor two moving from a far-southern location 
to a far-northern one, where it reaches almost 12%. If one asks for a minimal neutrino 
path-length of L = 3000 km 20^, the shadowing probability is about 10% less than 
the shown Earth-crossing case. However, the minimum path-length to detect Earth 
signatures depends on the detailed features of the flavor- dependent neutrino fluxes and 
on the detector properties. Therefore, to keep our discussion simple and general we 
restrict ourselves to the illustrative cases of Earth and core crossing. Note that the 
generic case can be calculated with the online tool |2E]- 

Some representative locations are tabulated in Table [T] The list is not exhaustive, 
but rather meant to provide representative and complementary geographic locations. It 
includes two existing locations of high-statistics experiments sensitive to SN neutrinos, 
i.e. Super-Kamiokande in Japan j32j and IceCube that is under construction at the 
South Pole ^^, and a few possible locations of next-generation detectors proposed 
in the literature. In particular, the possibility of a Megatonne water- Cherenkov 
detector is discussed worldwide, including the Underground nucleon decay and Neutrino 
Observatory (UNO) with a possible location in the Soudan mine (Minnesota, USA) j37j, 
the Hyper-Kamiokande (HK) detector in the Tochibora Mine (Kamioka region, Japan) 
[SB], and the MEgaton class PHYSics (MEMPHYS) detector at the Frejus site (France- 
Italy) jSn]. Different locations are under study for the large- volume scintillator detector 
LENA. The Pyhasalmi mine in Finland and Hawaii are of interest for geo-neutrino 
studies. Hawaii and Australia have also been discussed as favored sites for detecting the 
cosmic diffuse SN neutrino background because of the low reactor background [lOj. 

3.2. Two detectors 

Next we consider the simultaneous SN neutrino detection by two different detectors. A 
single detector can observe Earth matter effects only if it has good energy resolution and, 
of course, enough statistics. A scintillator detector like the 50-kton LENA project has 
good energy resolution and can detect unambiguously these signatures, whereas water 
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Cherenkov detectors have an intrinsically poorer energy resolution so that Megaton- 
scale masses would be necessary j^U]- However, smaller detectors, even if they can not 
resolve directly the modulations, could still reveal Earth effects by the comparison of 
the signals from a shadowed with an unshadowed detector. Of particular interest are 
Super-Kamiokande in Japan and IceCube at the South Pole that will be completed 
in a few years. In this case the Earth effect will be detected as a difference in the 
signal normalization between the two detectors -BJ. To this end it would be crucial that 
one of the detectors is shadowed whereas the other is not, i.e. we are interested in the 
probability for exactly one of them to be shadowed. 

The special location of IceCube makes this kind of problem a very simple 
generalization of the previous one-detector case. The shadowing condition of Eq. (fT3jl 
is satisfied or missed at the South Pole with a probability given by the lower or upper 
sign of Eq. ()19j) . and at any other location according to Eq. (j^^ or its complement to 1. 
Note that we can safely use the formulae of Sec. 13.11 since the location of one detector 
at the South Pole implies that no longitude difference enters the problem. 

In the general case of two detectors at geographic latitudes A and A' with a difference 
in longitude of A, the altitudes a and a' of an object located at (a, 6) are 

J sin a = sin A sin (5 + cos A cos (5 cos iJ , , , 

I sin a' = sin A' sin 5 + cos A' cos 6 cos{H + A) . 

To maximize the chance to detect Earth matter effects, the most interesting case for two 
general detectors is that exactly one of them sees a SN from a shadowed position, thus 
allowing one to compare a shadowed signal with an unshadowed one. The probability for 
detector I to be shadowed and detector II not is: (i) zero if one of the two probabilities 
vanishes when calculated as described in the previous section; (ii) trivially equal to 
the one-detector case when the probability of one of the single-detector conditions is 1; 
(iii) in all remaining cases it can be evaluated from the following system of inequalities 
where we omit (mod 2n) 

J — aiccos hfi{X, 6) < H < arccos/iK(A, 5) , , , 

I arccos h^iX', S) < H + A < 27r — arccos h^{\' , 5) . 

The fraction of the 27r-angle satisfying the previous conditions is the desired probability. 
The complementary situation, i.e. detector II is shadowed and detector I not, is 

J arccos /Ik(A, 5) < H <2ti — arccos h^{\, 5) , , , 

[ — arccos hi^{\' , 5) < H + A < arccos h^{X', 5) , 

apart for the trivial cases (i) and (ii) described above. 

It is also worthwhile to consider the situation that two detectors with excellent 
energy resolution are available, e.g. two scintillator detectors. We have already 
mentioned that next-generation scintillator detectors may be built at two different sites, 
e.g. one in Europe and another in Hawaii, for the purpose of geo-neutrino research. In 
this case, the most interesting case is that at least one detector sees a shadowed signal. 
The generalization of the previous formulae to this case is straightforward. 
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Our results are reported in Table |21 for the Earth-crossing case and in Table 01 for 
core crossing. The quantities listed in each cell at row i and column j are the following 
probabilities: 

P{i,j) Detector i shadowed, j not shadowed ; 
P{i,i) Detector i not shadowed, j shadowed ; 
P{i,j) Both detectors shadowed . (30) 

For example, in Table |21 in the row "Pyhasalmi" and column "South Pole" we read that 
the probability that a SN will be shadowed at Pyhasalmi but not at the South Pole 
is 0.519, for the reverse case is 0.353, whereas the chance that both detectors see it 
shadowed is 0.062. 

For the special case where one of the sites is the South Pole, it is intuitively obvious 
that Pyhasalmi is the best location among the proposed ones, being farthest at north. 
For the most conservative scenario, in which no new large detectors will be built, the 
combination of IceCube at the South Pole plus Super-Kamiokande detector already 
existing in Japan offers a 73% probability that a comparison between a shadowed and an 
un-shadowed neutrino signal could be observed, and 17% probability for this comparison 
in the case of core shadowing. Notice that the case of a detector in Pyhasalmi and 



Table 2. Shadowing probability for two detectors. The quantities reported in each 
cell are explained in Eq. piljl . 



LOCATIONS 


South Pole 


Sydney 


Hawaii 


Kamioka 


Frejus 


Soudan 




.519 


.457 


.285 


.179 


.065 


.148 


Pyhasalmi 


.353 


.315 


.231 


.157 


.052 


.139 




.062 


.130 


.296 


.400 


.516 


.433 




.475 


.437 


.187 


.221 


.162 




Soudan 


.317 


.310 


.142 


.208 


.158 






.097 


.135 


.386 


.351 


.410 






.461 


.484 


.326 


.230 






Frejus 


.307 


.361 


.285 


.220 








.107 


.084 


.242 


.339 








.435 


.278 


.184 






Kamioka 


.290 

.124 


.164 
.281 


.152 
.375 










.362 


.251 






Hawaii 


.249 
.165 


.168 

.277 












.160 






Sydney 


.129 

.285 
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Table 3. Same as Table El for core shadowing. 



LOCATIONS 


South Pole 


Sydney 


Hawaii 


Kamioka 


Frejus 


Soudan 




.115 


.116 


.115 


.116 


.041 


.110 


Pyhasalmi 


.064 


.069 


.082 


.114 


.035 


.107 




.000 


.000 


.000 


.000 


.074 


.005 




.112 


.112 


.107 


.112 


.111 




Soudan 


.064 


.069 


.077 


.104 


.108 






.000 


.000 


.005 


.000 


.001 






.110 


.110 


.110 


.110 






Frejus 


.064 


.007 


.082 


.104 








.000 


.000 


.000 


.000 








.104 


.104 


.100 






Kamioka 


.064 
.000 


.069 
.000 


.078 
.003 










.082 


.082 






Hawaii 


.064 
.000 


.069 
.000 












.062 






Sydney 


.057 
.006 













another one in Hawaii, which is of importance for geo-neutrino research, offers also a 
nice opportunity for Earth effect detection in SN neutrinos. The probability that exactly 
one of them is in a shadowed position exceeds 50% (20% for core shadowing), whereas 
the probability that one or both are shadowed is about 80% (20% for core shadowing). 

4. Conclusions 



The possibility to detect Earth matter effects in the neutrino signal from the next 
Galactic SN is a powerful tool to probe the neutrino mass hierarchy. Next-generation 
large volume detectors with excellent energy resolution would offer the opportunity to 
detect directly the specific signature of an energy- dependent modulation of the measured 
neutrino flux. Even if current detectors could not reconstruct directly these modulations, 
the comparison of the neutrino signal in a shadowed with an unshadowed detector could 
allow one to detect Earth effects. Of course, it is assumed that the location of the SN 
in the sky can be determined by observations in the electromagnetic spectrum. In the 
unlikely case that this is not possible, neutrinos alone can be enough to determine the 
SN location with sufficient precision |¥T1I321- 

Motivated by these opportunities, we have provided the ffist detailed study of the 
probability that a detector at a given geographic latitude will observe a Galactic SN 
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in an Earth-shadowed or core-shadowed position. We have shown that this probabihty 
is rather insensitive to detailed assumptions about the uncertain distribution of core- 
collapse events in the Galaxy. The main effect is the simple geometric constraint that 
SNe occur within the Galactic disk. 

We find that a far-northern location such as the Pyhasalmi mine in Finland, the 
preferred site for the LENA scintillator detector, is almost optimal for observing a SN 
signal shadowed by the Earth. The shadowing probability is close to 60%, against an 
average of 50% for a random location on the Earth. The shadowing probability Po{X) 
depends only mildly on the latitude A for A > 40°, a condition fulfilled by most of the 
northern locations proposed for next-generation experiments. The behavior for core 
shadowing is similar, although the advantage of a northern site is more pronounced. 
The core-crossing probability varies by almost a factor of two between a far southern 
and a far northern site, where it reaches almost 12%. 

We have also studied the case of two detectors. It is of interest because one may be 
able to compare a shadowed with an unshadowed SN neutrino signal and thus diagnose 
Earth effects even if both detectors lack the energy resolution that is necessary to observe 
a modulated signal. On the other hand, if both detectors see a shadowed signal one 
may perform Earth tomography because the observed neutrinos would cross different 
geophysical layers jlHl- For the South Pole, where IceCube will be completed in a few 
years, combined with Kamioka the probability that at least one of them is shadowed is 
almost 75%, for the South Pole combined with Pyhasalmi it is almost 90%. 

One particular scenario consists of a large-volume scintillator detector located at 
Pyhasalmi to measure the geo-neutrino flux in a continental location and another such 
detector in Hawaii to measure it in an oceanic location. The probability that exactly 
one of them is shadowed exceeds 50% whereas the probability that one or both are 
shadowed is about 80%. Therefore, Pyhasalmi and Hawaii are not only complementary 
for the purpose of geo-neutrino observations, but also for observing Earth matter effects 
in SN neutrinos. 

The various probabilities relevant for two detectors at arbitrary locations are 
difficult to represent in a useful figure or table. Therefore, based on our benchmark 
Galactic SN distribution we provide an online tool that allows one to calculate the 
various Earth and core-shadowing probabilities for one or two detectors at arbitrary 
geographic locations ^Hj . 
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